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Abstract 



In this paper we study the large-A'^ Umits of the integrable A^-state chiral Potts 
model. Three chiral solutions of the star-triangle equations are derived, with states 
taken from all integers, or from a finite or infinite real interval. These solutions 
are expected to be chiral-field lattice deformations of parafermionic conformal field 
CSJ . theories. A new two-sided hypergeometric identity is derived as a corollary. 
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^ . 1 Introduction 

M'. 

When the integrable A^-state chiral Potts model was introduced, it was the first 
example of an exactly solvable lattice model whose Boltzmann weights both 
require the use of higher-genus algebraic functions for their uniformization and 
do not have "the difference property" [1-6]. Since then, much has been written 
about many aspects of this model and we refer the reader to the recent review 
[7] for more information. In this paper we shall concentrate our attention on 
just one aspect, namely the large A^ limit. We have written about this once 
before [8] , but we can now present a much more complete and improved version 
containing several new results in addition. 
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The chiral Potts model is a spin model on a two-dimensional lattice (or more 
generally a planar graph). At each site (or vertex) of the lattice (or graph), 
there is a state variable or "spin" that takes on N values a,b,- ■ ■ = 1,2,- ■ ■ ,N, 
(mod A^). The Boltzmann weights are associated with pair interactions along 
edges. We assume that there are two types of such weights W and W, which 
on a square lattice would correspond to horizontal and vertical interactions. 
We assume also that the weights only depend on the difference modulo A^ of 
the two spin states a and b at the two endpoints of each edge, which is the 
Potts property. The chiral characteiQ (handedness or breakdown of parity) is 
expressed by W{a — b) ^W{b — a) and can only occur if A^ > 2. 

The integrable chiral Potts model [1-6] is a nontrivial generalization of the 
critical Fateev-Zamolodchikov model [9]. The fact that its rapidity variables 
lie on a higher-genus curve makes this model special among the many solvable 
lattice models. In spite of this several results exist for it. Therefore, its large- A^ 
limits should provide interesting generalizations of certain nonchiral cxo-state 
models of Fateev and Zamolodchikov [9-11], very different from the SOS- 
model of Baxter [12] and the few other cxo-state models [13-16] that have 
been introduced. From existing thermodynamic results for the finite- A^ case, 
we can infer corresponding results for the N = oo cases that may be of interest 
in later studies. We expect, for example, a direct relation with new integrable 
chiral-field deformations of parafermionic conformal field theories. 

This paper is organized as follows. In section 2 we present the Boltzmann 
weights of the integrable chiral Potts model and its dual model, also adding 
new details not given in [8]. In section 3 we give the three different large- A^ 
limits of the weights, while treating the more technical details in Appendix 
B. The three corresponding large- A^ limits of the star-triangle equations are 
given in detail in section 4. In section 5, it is shown that the results of the 
previous section 4 imply a new two-sided hypergeometric summation formula. 
Finally, a short discussion is given in section 6. 



2 Integrable iV-state Chiral Potts Model 



In this section we shall review earlier results on our higher-genus solution of 
the star-triangle equations for the chiral Potts model [1-6] and present in more 
detail a reparametrization [8,17] that is particularly suitable for the large- A^ 
limit. 



^ This chiral aspect allows us to mimic the effect of further-neighbor interactions 
within the context of a nearest-neighbor interaction model, see [7] and references 
quoted there. 
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Fig. 1. The square lattice represented by solid lines is diagonally oriented here, 
with spins taking N different values at its vertices and Boltzmann weights W and 
W associated with pair interactions along the two types of edges. The positions 
of the spins of the dual lattice are indicated by open circles. The medial graph is 
represented by the dashed oriented horizontal and vertical lines, also called "rapidity 
lines," carrying the spectral or rapidity variables p and q. 

2. 1 Star- Triangle Equation for Chiral Potts Model 



The iV-state chiral Potts model can be defined on a general graph with spin 
states a,b, ■ ■ ■ taking values 1, ■ ■ ■ , iV on the vertices and Boltzmann weights 
W{a, h) =W{a — h) associated with edges. W{n) is periodic in n mod N . 

In the integrable model, one assumes that there are oriented straight lines (the 
rapidity lines) on the medial graph, which are dashed lines shown in Fig. 1 for 
the case of a square lattice. They are obtained by connecting the middles of 
all pairs of edges (solid lines in the figure) that are incident to a single site and 
share a common face. No more than two rapidity lines meet at any given point. 
These lines carry variables p, q, ■ ■ ■ and arrows specifying their orientations. 
In nearly all solvable models the weights depend on the differences of these 
rapidity variables. For our class of integrable spin-pair interaction models the 
weights can be graphically represented as in Fig. 2. These weights must satisfy 
the star-triangle equation 
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i?pqr W^pq(a - 6) W ^r{h " c) W^r{a - c). 



(2.1) 
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Fig. 2. Boltzmann weights Wpq(a — b) and Wpq{a — b). We need to put an arrow on 
each edge to distinguish Wpq{a — b) from Wpq{b — a). Note the relative orientation 



of this arrow with respect to the orientations of the two rapidity hues in each case. 
Here the factor R^qj. can be determined as [5,6,18] 
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(2.2) 
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_ 2TTi/N _ 2-K^/^/N 



(2.3) 



The easiest way to derive (2.2) is to set a = in (2.1) and then to take the 
determinant with respect to the matrix indices h and c, leading to determinants 
of products of diagonal and cyclic matrices; this argument first appeared in 
print in [18]. The star-triangle equation (2.1) can be symbolically represented 
as in Fig. 3. 
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Fig. 3. The star-triangle relations, which allow one to move a rapidity line p through 
a vertex, which is the intersection of two other rapidity lines q and r. 



2.2 Weights of Integrable Chiral Potts Model 



In [1-6] one family of integrable chiral Potts models, with an arbitrary number 
of states per site N > 2, has been deduced with weights Wpq{a — h) and 
W^pqlct — b) satisfying the star-triangle equation (2.1) for all a, 6, c = 1, ... , A^. 
These weights are given by 



W^pq(O) \iij fj^yp-x^ui' 

W^pq(O) ^'^"'^"^ /i Z/q-l/p^^ ^ ' 

Here,[^the parameters p = (xp, ?/p, /ip) and q= (xq, |/q, yUq) are restricted by the 
two periodicity requirements PVpq(A^ + n) = W-pc^{n), iypq(A^ + n) = W 
yielding 



pq 



n] 



/^p)^_ 


yN _ N 

yp "f-q 

yq 3;p 


(/^p/^q)^ = 


< - < 
^p ^q 



(2.5) 



which can be recombined as 

l±^^^ 1±< 



\±, 



(2.6) 



independent of p and q. We write \± = ±c{l^k')/k, with c a constant that can 
be absorbed by a trivial rescaling of all the Xp,yp,Xq, yq by a common factor 
that will drop out of (2.4), and with k and k' numbers related by k +k' = 1. 
Then the conditions (2.5) reduce to 

f,^ = k'/{l-kx^) = {l-ky^)/k', < + yp^ = A:(l + «)- (2.7) 



These equations describe a complex curve, which is the intersection of two 
"Fermat cylinders," and the genus of this curve is ^f = N'^{N — 2) + 1. For each 
Boltzmann weight the two line (or rapidity) variables p and q are two points on 
this higher-genus algebraic curve, so that the usual difference-variable trans- 
formation cannot be carried out, except for special subcases where the genus 
degenerates to g < 1. Here the substitutions k = 0, k' = ±1 reduce the curve 



^ To connect with the original homogeneous notation [5,6], we must set Xp = Op/dp, 
ijp = bp/cp, /ip = dp/cp, and similarly with p replaced by q, r, • • •. A proof that the 
star-triangle equation (2.1) is satisfied is given in the appendix of [6]. 



(2.7) to a product of genus- zero curves and the weights degenerate to those of 
the self-dual Fateev-Zamolodchikov model [9]. 

In general the model is not self-dual and the dual weights are found by Fourier 
transform. We first note that the periodic weights in (2.4) are of the form 



so that the linear recursion relation 

(X4 - X3cu")iy(n) = (xi - X2u'')W{n - 1) (2.9) 

is satisfied. We can apply the Fourier (or duality) transformation^ 

N~l 

W^^\m) = N-^ Y. ^""""W^H, (2.10) 

n=0 



leading to 



or 



= uj-"'xiW^^\m) - uj^-"'x2W^^\m - 1), (2.11) 



ly (f)(0) f}^ xi-x^uj^ ■ ^ ' ' 

Therefore, the weights dual to (2.4) are 

W'il{0) j=i ^/^p/iqa;p - VqUJ^ ' 
which are again both of the form (2.8). 



^ Here we added a normalization factor A^ ^ that is needed in the A^ — > oo limit. 



2.3 Reparametrization 



In order to proceed, we introduce new parameters to describe the higher-genus 
curve of rapidity variables. These parameters are real when the Boltzmann 
weights Wpq{a — h) and iypq(a — h) are real and positive. We begin with the 
substitutions^ 



a;„ = e-p'-, 2/p = cu^e'^'=/^, Xq = e''^'^/^, i/q = c^5 e'^A^, (2.14) 






SO that from the last identity in (2.7) we find 



^^ ^ sin i(gp - 0p) ^ sin \{e^ - 0q) 
sini(^p + 0p) sini(^q + 0q)' 



This is equivalent to 

e'^p +k 
l + k e'^p 



-''' = ^^^' (2-16) 



and similarly with p replaced by q. From (2.16) we have 

2fc + (l + fc^)cosgp . (l-fc^)singp 

cos0p= . .,2,nu a~' sm0p= ^ ^.^ . ^, ^- 2.17 

1 + K^ + 2A; cos t^n 1 + K^ + 2A;cos6'p 



We will also need two parameters that will describe the dual model, see e.j 
(2.29). The first one is given by 



6'p + 0p 1 sin 6'p 

= — arctan — - 

2n TT cos 9^ + k' 



Ap = ^JL^ = ± arctan -£^P^, (2.18) 



where the last step follows from (2.16). Also using (2.15) we find 

6'p — 0p = 2 arcsin(/c sin vr Ap) , (2-19) 

6'p = vrAp + arcsin(A;sin7rAp), 0p = vrAp — arcsin(/csin7rAp), (2.20) 

which expresses 6'p and 0p in terms of Ap. The other parameter jp is defined 

by 

piiSp it. 
g7p±^iAp ^ (2.21) 



^ Our definitions of 6p and (/>p differ by a factor A^ from Baxter's [17]. This change 
of normalization will be necessary in the large A^ limit. 



These two expressions are equivalent in view of (2.16) and (2.18). Multiplying 
them and using the second equality in (2.17) we find 



g27p _ l + A:^ + 2A;cosgp ^ singp ^^ 22) 

I — k^ sin ( 



'p 



Because of (2.20) this 7p is also a function of Ap, i.e. 



,/ 1 — A;2 sin^TrAp ± k cos vrAp 
:,±7p — J 1 

From (2.7) and (2.14) we have /i^ = (1 + A; e^p)/A;' so that 

^ e"^p sm 9 ^ 



(2.23) 



^' = Uri^j • P-2*' 



after using (2.18) and (2.21). 

With the help of (2.14) and (2.24) we can now rewrite the results (2.4) as [17] 



W,^in) ^ / singpsin0q \ "^'^ " sin[7r(j - i)/iV - (g^ - 0p)/2iV] 
W^pq(O) I^sin0qsin0p j /J sin[7r(j - i)/iV + (0^ - ^p)/2iV] ' 

(2.25) 

^pq(n) ^ / singpsingq y/'^ Jt sin[7r(j-l)/iV+(c/,q-./,p)/2iV] 
Wpq(O) Um0psin0j /J sm[7Tj/N-{e^-e^)/2N] 

(2.26) 
Similarly, using (2.14) and (2.24), their Fourier transforms (2.10) become[^ 



^ffH ^ ^in(<Ap-.p+^,-.,)/2;v A sin[7r(j-l)/iV+(0q-0p)/2iV] 
W^?(0) /=i sin[7rj7iV-(0~q-0~p)/2iV] ' 

(2.27) 

^S(^) ^ ,in(.p-,p-.,+,,)/2;v A sin[7r(j-l)/iV-(0q-gp)/2iV] 



e 



n 



M^S(O) ~' /ii sin[7r(j - i)/iV + (^^ - 0p)/2iV] ' 

(2.28) 



^ Eq. (7) and (8) of [8] have misprints, which can be corrected by replacing n by 
A^ — n in their left-hand sides. 



where 



0p = TrAp - i7p = K^p + 0p) - |i log — — ^, 

_ sin ^ 

^p = TrAp + i7p = K^p + 0p) + \i log --f . (2.29) 

By direct substitution we can show that if the weights satisfy the star-triangle 
equation (2.1) then their Fourier transforms satisfy the star-triangle equation 



^ Wma)Wi^Xh)W^l{^ + h) 



■ • qr v-y ■ ■ pr \" / pq 
-'I'pqr 



Af-1 

E 

d=0 



Y. W^S(& - rf) W^£H« + h-d) W^Xd)- (2.30) 



This equation has the exact same form as equation (2.1), as can be seen 
replacing a-^a — b,b^b — c, a + b^a — c, and c + d ^ d. Therefore, from 
the proof [6] that the weights (2.25), (2.26) satisfy (2.1) we conclude that the 
weights (2.27), (2.28) satisfy (2.30). 



For 6'p = 0p, 6'q 


= 0q we 


recover 1 


he 


solution with 








W^«(n) 


W{n) 


W{N- 


n) 


w^m(o) 


W{0) 


W{0) 




W^'\n) 


W{n) 


W{N- 


n) 



w^m(o) w{<d) w{<d) ' ^^'^^^ 

which are trigonometric expressions {g = 0) of the difference variable 6q — 6p. 
In this nonchiral special case the Boltzmann weights depend only on this one 
parameter, which is the difference of two rapidity variables. The more general 
chiral weights depend on the two rapidity variables separately, living on a 
higher-genus curve. 



3 The A^ ^ oo Limit of the Boltzmann Weights 



In this section we shall obtain the N —^ oo limit of the Boltzmann weights of 
the previous section. We shall give explicit formulae for all three regimes. 



3. 1 General form of the Boltzmann weights 



Note that the Boltzmann weights (2.25) and (2.26) or their dual weights (2.27) 
and (2.28) all have the product form 

W{n) _ .- sin{iT{j + a-l)/N) 

W{0) l{sm{n{j+(3-l)/Ny ^'^ 



where 

A = sin TT/S/ sin na, (3.2) 

with a and P given constants depending on parameters 9p, 9^, 0p, and 0q 
satisfying (2.15). Also, the condition on A guarantees that W{n + N) = W{n), 
using a trivial exercise on complex exponentials. 

More precisely, we have to use in case of (2.25) and (2.26) 



± yp t^q n ^ I yq '^p /o o\ 



0q 


-0, 




27r 


^ 


-^q 



PI 2 27r ' ^P'l 

_ rq ~ rp 75 i , "P "q ro a\ 

P'l " 27r ' /^Pq = 1 + ^^' ^3.4) 

which all four satisfy equations of the form ^pq + .^qj. — ,^pr = ^qq, or 



"pq + "qr - "pr = /^pq + /^qr " /^pr> (3-5) 

"pq + "qr = "pr, ^pq + ^qr = 1 + ^pr' (3-6) 

Similarly, we have to use in case of (2.27) and (2.28) 



«S = ^^ ^2 = ^ + ^"' ^'-'^ 



"pq 2 27r ' '^P^ 2 27r ' ^ 



1 Or, — (J)n —/r\ 1 Un- 
satisfying 

«S + «ff = <^ P^vl + /^S = 1 + /^?r\ (3.9) 
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The four corresponding constants A, A, A^^\ A^^\ as given in (3.2) with the 
corresponding a and [3 substituted, are worked out in Appendix A and they 
agree with (2.25) through (2.28), as was to be expected. 

Important symmetries of weight (3.1) are 

W{n\a,(3) W{0\(3,a) W{n±N\a,(3) 



W{0\a,p) W{n\p,a) W{0\a,p) ' 

W{-n\a,(3) W{N-n\a,(3) W{n\l - (3,1 - a) 



W{0\a,(3) W{0\a,f3) W{0\1 - /3,1 - a) 



(3.11) 
(3.12) 



which is easily verified from (3.1). This allows us to restrict ourselves to study 
W{n) only for < n < ^N, while reducing the other case ^N < n < N or 
equivalently —^N < n < to this case. This symmetry shows up explicitly in 
the following, particularly in (3.17), (3.19), (3.21), and (B.4). 

We can conclude from (3.12) that W{N — n) = W{n) for a + j3 = 1. Then 
the chirality disappears and the model reduces to the model of Fateev and 
Zamolodchikov [9]. 



3.2 General iV — ^ oo formula 



Naively, in the limit A^ — > cxd, we can drop the sin symbols in (3.1). This leads 
us to introduce the function 



T{(3 + n)T{a) /= ^ J + /9 - 1 (/9)„ 

= n^^ = |l^^, ifn<0, (3.13) 

where r(a;) is the Gamma function and {x)n = T{x+n)/r{x) the Pochhammer 
symbol [19,20]. The finite- A^ corrections are described by the function 

^ , , , " sin[7r(j + a-l)/A^] ,^ _, ^, 

S^ia) = log n ^^T TTT^, 3.14 



which has an asymptotic expansion derived in Appendix B. Using (B.2) there, 
we immediately have an asymptotic expansion formula for (3.1) in terms of 



11 



powers of 1/A^, i.e.Q 



log 



W{n) 



log U"/^P(n|a,/3)l + Sn{a) - 5„(/5) 



1^(0) 

oo 

= log[A"/^P(n|a,/?)j +J2 



1=0 



X 



^ d \ , fsinz 



B;+i(a)-B;+i(/?) /7ry 
(/ + !)! ViV 



z=TTn/N 



i) '- 



smz 



2=0 



(3.15) 



where the Bm(a;) are Bernoulli polynomials [19,20]. Only the term I = will 
be relevant in the limit N ^ oo and the terms / > 1 are finite- A^ corrections, 
for which bounds are derived in Appendix B. Using Bi(a;) = x — | and B2(x) = 
x^ — X + I [19,20] and restricting ourselves to / < 1 we can rewrite (3.15) as 



W^H_^n/A.p(^|^^^/sin(WiV)' 



a-/3 



W{Q) 



X exp 



V TTTl/N 

2iV 



nn N 
cot — 



0(Ar-2) 



(3.16) 



The last line of (3.16) gives the leading correction for large N and can be 
ignored in the limit. We can use (3.16) to study three regimes for the large N 
limit. We shall work this out in the following three subsections. 



3.3 The regime I: N ^ oo, n finite 



First we study the limit A^ -^ oo, while n remains finite. In this case, (3.16) 
results in 



W{n) _ T{a + n)T{p) _ r(l - /3 - n)T{l - a) 
W{0) ~ T{(3 + n)T{a) ~ T{1 - a - n)T{l - (3)'' 



— oo < 77. < oo. 



(3.17) 



which is just the naive limit P{n\a,f3) given in (3.13). Using the well-known 
asymptotic expansion formula of the Gamma function [19,20], we have 



T{z + a) 

r(^ + /?) 



^ (l + (^ - ^)('^^+ ^ - ^) + o(.-)) , (3.18) 



^ Formula (12) in [8] has the higher orders misprinted and is only correct to the 
order needed in the actual N ^ oo limits, which are presented correctly in [8]. 
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for z -^ +00, which is also equation 1.18(4) of [19]. Therefore, 

= H"'^ r[il^) (1 + O(l^r')) > for n ^ -00. (3.19) 

This shows that the Bohzmann weights vanish in the hmit whenever SfJo; < 
where 3fJz is the real part of z. 



3.4 The regime II: N,n ^> 00, n/N finite 



For the second regime we study N, n ^ 00 such that 

27m 
x^^ (3.20) 



remains finite. Consequently, the weights W{n) in (3.1), which originally took 
N different values and which were periodic modulo A^, now depend on the 
continuous spin values x and they are periodic modulo 27r. 

We can now substitute the asymptotic formula (3.19) for P{n\a, 13) into (3.16), 
while assuming without loss of generality ~\N < n < ^N and rearranging 
the resulting expression. We immediately arrive at 






— /A^ 1 \"~^ T(l- a) 
= W{0)A^-(-sml\x\) -)^ -J-, ii-l7i<x<0. (3.21) 

This can be summarized as a function periodic modulo 27r, i.e. 

W{x) = CA^^l^\ sinlx'"''^, (3.22) 



where \_x\ stands for the largest integer < x and 

C-W(0)(^Y''^^ sinvr^ _ rXa)ra^a) 



In this regime II, we have to rescale W{0) with a power of A^ as A^ ^ 00 in 
order to keep the constant C finite. 
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For the special case a + f3 = 1, we have A = 1 and the chirahty vanishes. This 
special limit has been mentioned first by Fateev and Zamolodchikov in [9] and 
is generalized above to a + (3 ^ 1. 



3.5 The regime III: N,n -^ oo, n/N -^ 



A crossover regime intermediate between regimes I and II appears when both 
iV, n — > oo such that n/ip{N) — > x for some function v'(A^) with v'(A^) — > oo 
and f{N)/N -^ 0. We have 

W{x) = D A'^"^°^"^ l^r"^, -cx)< X < oo, (3.24) 



which is a chiral generalization of the Boltzmann weight in Zamolodchikov's 
Fishnet Model [10]. Here, 

D = W{Q)^{Nr~^^Al, (3.25) 



implying again that W^(0) need be suitably rescaled in the limit A^ — ^ oo. We 
note that (3.24) is also the asymptotic large-n behavior (3.19) for regime I 
and the small-a; limiting behavior of (3.22) in regime II. The sign function 
in (3.24) arises as coefficients in (3.19) and (3.21) differ by a factor A for n 
positive or negative, see also (3.23). 

We note that we can reproduce the previously known cases [9,10] by setting 
a + {3 = 1 and A = 1. Now we have only one condition (3.2) on A, i.e. 
A = sinTT/S/sinTTQ;. This provides us with the deformations (3.17), (3.22), and 
(3.24), which define integrable field theories with chirality. 



3. 6 Duality of regime I and regime II 



The limiting Boltzmann weights in Regimes I and II are each other's dual 
under Fourier duality transformation. More precisely, if the limiting weights 
are in Regime II, their Fourier transforms are in Regime I, and vice versa. This 
follows from the way that we have constructed the limits. However, there is also 
a direct way to show this, as the infinite Fourier sum can be performed using a 
transformation formula of the Gauss hypergeometric function F(a, 6; c; x) [21]. 
Thus we obtain a formula for the double-sided hypergeometric function iHi 
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as defined for example by Slater [22], i.e.Q 



iHi 



- E 



a, 



(/3)n 



" inx 



iHi 



1-/5 
1 — a 



F(a, 1; /3; e^") + F(l - /5, 1; 1 - a; e"'") - 1 
2/3— ir(l-a)r(/?) ,n_._,K._. 



T{P-a) 



^.(i_,_^)(,_^)/2 (sinl^)/3— 1^ 



(3.26) 



for < X < 27r (and periodically extended mod 27r). The inverse Fourier 
transform of (3.26) corresponds to integral 3.892.1 of [20], where we need to 

replace ^/^— *>/? — a, /5^— *>« + /? — 1 + 2n, x \-^ n — ^z, resulting in 



271 



2l-/3+a gi7r(l-a-/3)/2 f_l\n 

2-kJ ~" V— -2-J ~ (/5 _ a)B(/5 + n, 1 - a - n) 



dz e-'"^+'(i-"-'3)^/2 ^gi^ 1^) 



1 ./3-Q-l 



r(/? + n)r(l - a - n) ' 



(3.27) 



in agreement with (3.26) as 
, , r(x + n) 



rri-x) 



-IT 



[X 



T{x) 



r(l — a; — n) (1 — x)_„ 



(3.28) 



It is easily verified that, up to a possible overall constant factor, (3.26) and 
(3.27) relate the regime I result (3.17) with the regime II result (3.22), provided 
an appropriate transformation of the a and (3 parameters is made. From (3.3), 
(A.l), (2.29), and (A. 3), we find for this duality transformation 



/3-a-l =a^^^ -I3^^\ 
a + (3-1 = --logA(f), 



^(f) 



smvr 



/?(f) 



71 ' sinvra^^)' 

whereas, from (3.3), (A.l), (2.29), and (3.7), its inverse is given by 



(3.29) 



a«+/5«-l = -logA 



A 



sinvr/^ 



smvra 



(3.30) 



s Start from 2.9 (27) of [19], with a = a, b = 1, c = (3, z = e^^, substituting the 
definitions 2.9 (1), (13), (22). Next use eqs. 2.1.2 (6) and F{a,0;c;z) = 1. 
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The equations (3.29) and (3.30) differ only by one minus sign in front of the i, 
coming from the corresponding difference of the finite Fourier transform and 
its inverse. 



4 The A^ ^ oo Limit of the Star-Triangle Equation 



In this section we shall examine the various limits of the star-triangle equation 
(2.1), now we have obtained explicit prescriptions on how to take the N —^ oo 
limit of the Boltzmann weights of the A^-state chiral Potts model. 



4-1 Principal domain 

There are several parameter domains that require separate treatment. But 
from now on, we shall assume that all W{n) encountered become vanishingly 
small whenever n,N — n ^ 0. In view of (3.19) and (3.21), this means that 
all 3fJ(/3 - a) > 0. 

From (3.3), (3.4), (3.7), (3.8), (3.29), and (3.30) we find 

f3-a = l-^ + a = l- p^^^ + a(^) = 'p^^^ - a^^\ (4.1) 



independent of the choice of rapidity variables p, q, which we have suppressed 
in (4.1). Therefore, we define our principal domain by the condition 

< 3f?(/5 - a) < 1, (4.2) 



for all occurring (a,/3) pairs. This condition is easily satisfied, even with ttp^, 
/5pq, ctpj,, /3pj,, ciqj,, /3qj,, and their barred versions, all being real. 

The summation over d in (2.1) has to be split in several pieces as we must 
choose to which of the three N ^ oo regimes each of the three weights in the 
summand belongs. We shall see that under condition (4.2) all pieces but one 
can be ignored and that the three types of large-A^ behavior I, II, or III do 
not mix: If we take the three spin states a, b, and c in (2.1) such that all three 
weights in the right-hand side of (2.1) are in the same regime, the dominant 
part of the sum over spin state d comes from the piece with all three weights 
in the left-hand side of (2.1) being in the identical regime. 
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4-2 The constant R 



pqr 



The next calculation to be done is the evaluation of the large- iV limit of the 
constants -Rpqr or Fpq given in (2.2). We break this up in several steps. 

First, we can rewrite the Fpq in (2.2) as 



^m = N^^eMm-Lr>.), 



W^pq(O) 



(4.3) 



where 



1 



N 



w^pq(0 



^pq-^Elog-^^^(Oy 



f(f) 



N 






(4.4) 



Here both the Lpq and the L^). can be evaluated in an identical fashion under 
condition (4.2), substituting the regime II asymptotic form (3.22) and (3.23), 
while replacing the sum by an integral. Therefore, in an obvious simplification 
of notation suppressing the rapidity subscripts, Lpq becomes 



L = log[{N/nr'^r{f3)/r{a) 

1 N 

+— ^ [[(n/N) - [n/N\) log A + (a - P) log sin |7™/A^| 
^\og[iN/nr-^TiP)/Tia) 



2n 



2tt 



dx 



[xJI-k) - La:/27rj) log A + (a - /?) log sin \\x\ 



\oA{Nl2Tiy<' A^I''Y{i5)lY{a] 



for iV — ^ cxD, 



(4.5) 



where the elementary integral 4.224.3 of [20] has been used. Similarly, 

m ^ log[(iV/27r)"**'-'^'' A(f)'/' r(:^m)/r(a^'))l , for iV ^ oo. (4.6) 



It is easily checked that the corrections to (4.5) and (4.6) are irrelevant in the 
large- A^ limit. 

Similarly, W^*^^)(0)/W^(0) is also dominated by the regime II contribution given 
by (3.22) and (3.23), i.e. 
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N _ _ _ 

Y^[{N/^f-^T{$)/T{a)\A^ 



(n/N)-\n/N\ i 



sinyrn/A^I 



a-l3 



'^-Prm\ 



{N/nr-'T{(3)/T{a 



1 



27r 



P/2^r.in i^^"-^ 



2tx J ^ 





We have used the integral (3.27) for n = 0, after substituting the barred 
version of (3.30), i.e. 



a 



f3 = /^m - cjW - 1, log A = TTi (1 - a(^) - /^W). 



(4i 



Then the last line of (4.7) is obtained using (4.8) once more. As before in 
(4.5) and (4.6), (4.7) gives the coefficient of the leading iV-power as A^ ^ oo. 
Correction terms can be obtained, but they will not be needed. 

We can now substitute (4.4)-(4.8) into (4.3) and simplify the result. This is 
worked out in Appendix C and the result is 



pq 



{N/2ir\ 



/3pq"apq .-1/2- 



^pq(o)r(apq)r(/3pjr(i-apq) 

W'pq(O) r(/3pq)r(^pq-apq) ' 



(4.9) 



giving us the desired expression for _Rpqr = -Fpqi^qr/-^pr- We note that the factor 
in square brackets in (4.9) cancels out in view of (3.5) and (A.l). Hence, 



lim i?, 



pqr 



iypq(0)_^pr(0)iyqr(0) 
W^pq(0)TFpr(0)l^qr(0) 



pqr' 



(4.10) 



where 



oo _ /pq/qr 

pqr 



/, 



/pq — 



r(«pq)r(/3pjr(i 



a 



pq; 



' Jpq 



pr 



r(/3pq)r(/3pq - apq) 



(4.11) 



We see that -Rpqr is independent of A^ in leading order in the large- A^ limit. 
Finite- A^ corrections can be worked out but are not needed here, as we shall 
only consider the actual N ^ oo limit. 

Note that /pq and r^j, are invariant under cJpq h-> 1 — /^p^, p^^ i-^ 1 — a^^, 
whereas Fpq and r|^j, are invariant under a^^ \-^ 1 — P^^, P^^ h-i> 1 — a^^, but 
/pq is not. 



4-3 Regime I 



We can first consider the large- A^ limit of (2.1), while keeping a — b, a — c, and 
thus also b — c finite. Without loss of generality, we can then restrict ourselves 
to considering the case a, b, and c finite. The three Boltzmann weights in the 
right-hand side of (2.1) take the regime-I form (3.17). The sum over d in (2.1) 
needs to be split up in the limit. In one part \d\ remains finite but can become 
arbitrarily large and the three weights in the left-hand side of (2.1) also take 
the regime-I form. The summand of (2.1) then decays as \d\'^ for \d\ -^ oo. 
From (3.17), and using (3.6) and (4.1), we see that 



= ttqj. — Pqj, + Ppj. — ttpj, — i + ttpq — Ppq = —2, l^.iZJ 

SO that the sum over d converges as X) MP • 

In the part of the sum for which \d\/N does not tend to zero, the three weights 
in the left-hand side of (2.1) belong to regime II and now we can use (3.21) to 
show that the summand scales as A^** = A^"^. As the total sum has A^ terms, 
this contribution vanishes in the limit. 

The contribution of the crossover regime III connecting regimes I and II also 
vanishes. To show this in more detail, we can split the sum over rf in a piece 

-A^2/3 < c^ < JY2/3^ g^^^ g^ pjg(.g JY2/3 < ^ < JY _ JY2/3^ j^^ ^J^g g^g^ pj^^^g^ 

the summand is bounded by \d\~ -behavior as shown above and in (3.24) and 
(3.25); therefore, the error made by replacing the sum with —N'^l^ < d < N'^/^ 
by a sum — oo < d < oo vanishes as 0(A^~^/^). For the second piece we can 
use (3.21), (3.24) and (3.25) to show that the summand is 0{N'^'^/^) and with 
less than A^ terms of order 0(A^~*^/^) its contribution vanishes as 0(A^~^/^). 

To summarize, all six Boltzmann weights must take the regime-I form (3.17) 

(«„„) (tt,^n) 

WAn) = jff, W,,{n) = '-ff, (4.13) 

(Ppq)n (Ppq)n 

solving the star-triangle equation 



J2 W^,{b-d)W^,{a-d)Wp^{d-c) 

d=—oo 

= Cqr W^pq(a " &) W^,{b - c) W^,{a - c), (4.14) 

provided the a and P parameters satisfy (3.3), (3.4), and (2.16). In (4.13) and 
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(4.14) we have chosen the normahzation VFpq(O) = W^pq(O) = 1. 

Equation (4.14) is related to the Dougall-Ramanujan identity [23]. We shall 
return to this in the next section. 



4-4 Regime II 



We can next consider the large- A^ limit of (2.1), while keeping (a — b)/N, 
(a — c)/N, and {b — c)/N fixed and nonzero. The three Boltzmann weights 
in the right-hand side of (2.1) now take the regime-II form (3.21). Again, the 
sum over d in (2.1) needs to be split up in this limit. Now the dominant part 
is the one with all three weights in the left-hand side of (2.1) belonging to 
regime II. Only when d is close to a, b, or c, one of the weights can be of the 
form of regime I or III. Because of the principal domain condition (4.2) these 
contributions can be ignored in the large- iV limit and the sum can be replaced 
by an integral as is done twice in subsection 4.2. 

In the previous section we have seen that the summand in the left-hand side 
of (2.1) scales as A^~^ in the large- iV limit, when the three weights belong to 
regime II. The right-hand side of (2.1) now scales as A^** = A^~^, since 



n, L^pQ Pno ~r "nr Pnr: ' or P 



pq ''^pq pr r^pr ' "qr r^ qr 



a 



pq 



/5pq + /5pr-«pr^ ^ + "qr~/^qr = -1- (^-IS) 



Therefore, it is natural to multiply (2.1) by A^ and to replace A^ ^J2d by 
(27r)~^ / dx, so that the star-triangle equation becomes 



27r 
— / dwW^riy - w)Wpr{x - w)Wpq{w - Z) 




= R^^, W^^ix - y) W^,{y - z) W^,{x - z), (4.16) 

as also follows after suitable rescalings of the Boltzmann weights and -Rpqr, i.e. 



i?^,,^^im AT-i/^pqr (4.17) 

and the two equations (4.20) below. Equation (4.16) has the solution 



l^pq(x) = e(^-^^^(^ 



WUx) = e^'^'-'^^^^^ 



f ) 



sin |a; 



sin ix 



(4.18) 
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using the normalization Cpq = Cpq = 1 in (3.22), together with (3.3), (3.4), 
(2.18), (A.l), and (A.2), which imply 



A = e'^p"^" a -8 = A - A 

4q = e"^^'^ «pq - ^pq = Aq - Ap - 1. (4.19) 

Here 7p and Ap are related by (2.23). If Ap < Aq < Ar < 1+Ap all six Boltzmann 
weights in (4.16) are real and positive and the parameters are in the principal 
domain. 

According to (3.23), the condition Cpq = Cpq = 1 implies 

^^^ r(/3pq) Vvry r(/?pq) 

Substituting this in (4.17), while using (4.10) and (4.11), we arrive at 

OOO _ /pq/qr f ^ 1 £(^pq)£(«pj . _ £(«pq)£(i^apq) ., .. s 

Jpr ^ 1 (apq)i (Ppq) TCL (^pq - ttpq) 



where the it and A^ factors have been redistributed using (3.5) and (3.6). 

It can be shown that the Boltzmann weights obtained by dropping the integral 
part [a;/27rj in (4.18) also satisfy the same star-triangle equation (4.16). The 
resulting chiral solution can be viewed as the nonchiral Fateev-Zamolodchikov 
large- A^ solution [9] with a site-dependent gauge transformation. 

Since W^pq(a;) and Wpq{x) as given in (4.18) are now functions of x periodic 
modulo 2tt, their Fourier transforms 



27r 2tt 

W^^?(J) = ^ / da: e-'^-- W^pq(a;), W^Hj) = ^ J dx e"'^'^ Wpq(x), 



(4.22) 

are over all integer value j, ranging from — oo to oo. Substituting (4.22) into 
(4.16), we find that these Fourier transforms satisfy the star-triangle equation 



^ WS{a-b)Wil\b-c)W^'i{a-c) 



poo " qr V"" -y ■ ■ pr V" ^y ' ■ pq 
^pqr 



E Wi2{b-d)W^;X<^-d)W^\d-c), (4.23) 



d=—oo 
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in which the sum is over all integer values of d. 



4.5 Regime III 



As the final case, we can consider the large- A^ limit of (2.1), while keeping 
(a — b)/ip{N), (a — c)/ip{N), and {b — c)/ip{N) fixed at finite and nonzero 
values. The three Boltzmann weights in the right-hand side of (2.1) now take 
the regime-Ill form (3.24). Also in this case, the sum over d in (2.1) needs to 
be split up. Now the dominant part is the one with all three weights in the 
left-hand side of (2.1) belonging to regime III. Only when d is close to a, b, or c, 
one of these weights can be of the form of regime I, while the other two weights 
are of the form of regime III. Because of the principal domain condition (4.2) 
these regime-I contributions can be ignored in the large- A^ limit and the sum 
can be replaced by an integral as is done twice already in subsections 4.2 and 
4.4. Also, the three weights in the summand of (2.1) can only simultaneously 
take the form (3.21) of regime II. This contribution is 0(A^~^) as in subsection 
4.3 and can be ignored. 

Hence, we can prove a star-triangle equation of the form (4.16), but with 
integration over (—00, 



+00 

J dwWqriy - W) Wpr{x - W) Wpq{w - z) 

—00 

= <r W^pq(a^ - y) W^y - Z) W^,{X - Z), (4.24) 

where 



^-, = lim ^{N)-'R^^,. (4.25) 



N^oo 



Equation (4.24) has the solution 



l^pq(x)=e-^^^''-^<')^'^"("V|'--"'' 



TFpq(x) = e -5(7P+7.)sign(.) |^|A,-Ap-1^ ^4_26) 

using the normalization Dp^ = D^q = 1 in (3.24), together with (4.19). Again, 
7p and Ap are related by (2.23). If Ap < Aq < Aj. < 1 -|- Ap all six Boltzmann 
weights in (4.24) are real and positive. 
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According to (3.25), the condition Dp^ = Dpq = 1 implies 

W^pq(O) = v^(iV)^--°- ^^^, TFpq(O) = v9(iV)^--°- ^i^. (4.27) 

r(ppq) r(ppq) 



Substituting this in (4.25), while using (4.10) and (4.11), we arrive at 

pco ^pq^qr j, __ Vf^pq/ V pqy n V pq/ V pq^ /a cyQ\ 

-'^pqr ~ p ? ipq — -^, N-r^/^ N •'Pq ~ ^,^ , ' l^-^°J 



where the </^(A^) factors have been redistributed using (3.5) and (3.6). 

It can be shown that the nonchiral Boltzmann weights obtained by setting 
7p = 7q = in (4.26) satisfy the same star-triangle equation (4.24). The 
resulting solution can be viewed as the Fishnet Model of Zamolodchikov [10]. 
It also generalizes Symanzik's conformal integral [11], which has been used 
by Zamolodchikov to prove the star-triangle equation for the Fishnet Model 
and which has also provided the proof for the N = oo Fateev-Zamolodchikov 
model [9] via a conformal transformation, C, = tan |-U7. 



4-6 Remark on R-matrix 



Following our joint work with Baxter [5], we can make an /?-matrix by taking 
the product of four weights in any of the above regimes I, II, or III. Taking 
four weights of type I, we have 

/?(a, 6, c, d) = Wp,q, (a - c) W^p.q, (c - b) Wp,q, {d - b) W^p,q, {a-d). (4.29) 



Similarly, we can also take four weights of type II or III. Then any such 
infinite-dimensional /?-matrix satisfies the usual Yang-Baxter equation. But 
these solutions are very different from those of [12-16]. 



5 Two-Sided Hypergeometric Sum 



In this section we shall rewrite the star-triangle equation (4.14) with solution 
(4.13) as a new double-sided hypergeometric identity. 
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5.1 More symmetric star-triangle equation 



The star-triangle equation (4.14) can be written in a more symmetric form 
with permutation symmetry among the Boltzmann weights in each of the two 
sides. 

We start by applying (3.12) to the W^^ in (4.14), followed by substituting 
d = —n and and replacing all six weights by (4.13). We arrive at 



J^ ("nr)^, ("r.r) , ( 1 " /^r-n ) , (^nn) u{Cir.r)u ("nr) 

V^ ^ qr^b+n ^ V^'a+n ^ '^V'i'c+n oo ^ Pq^g-fe ^ Pr^fe-c ^ 1^' a-c 

n=-oo (/?qr)b_^„ (/^pOa+n (^ ~ '^Pq)c+n if^P'Va-b (/^Pr)b-c {P^^)a-c 

(5.1) 

This result seems to invite us to introduce a more symmetric notation. For 
the three external spin states in the star-triangle equation we write 

mi = a, m2 = b, m^ = c. (5-2) 



It is logical to associate pr with 1, qr with 2, and pq with 3. The other 
quantities in (5.1) are then rewritten as[^ 



X\ — ^pr' 2 — ^qr' 3 — ^ '^^PQ' 

1/1 = /5pr) Z/2 = /5qr5 VS = '^ - "pq, 

Ui = Opr, U2 = ftqj., M3 = Opq, 

Vi = ^p„ V2 = /3qr, V3 = /3pq, (5.3) 

/l = /pr, /2 = /qr, /s = /pq, (5-4) 

where the fj for j = 1, 2, 3 are given by (4.11), i.e. 



„ r(xi)r(t;i)r(i-Mi) 

h - 



T{yi)T{vi -Ml) 



^^ T{u2)T{y2)T{l-X2) ^ ^^_ T{us)T{ysm-xs) ^ 
r{v2)T{y2 - X2) ' ^ r{v3)r{y3 - X3) ' 

and r^j, = J2f-il f\- Therefore, (5.1) takes the much more symmetric form 



^ The parameters Xj and y^, for j = 1, 2, 3, should not be confused with the rapidity 
variables Xp and y^ of (2.4). 
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E 



V-^l/mi+n V'^2Jm2+n V'^sJma+n 
wl/mi+n v"2J,n2+n wS/nia+n 

_ /2/3 (Mi)„j_„2 (m2)„2_„3 ("3)m,. 



■ni3 



^1 V 1/mi— m2 ^ ^ )m2^m'j, V 3/,ni- 



(5.6) 



"13 



The quantities (5.3) are not independent. From (3.3) and (3.4) we find that 
the six variables Uj and Vj depend in a hnear and symmetric fashion on the 
six variables Xj and Hj, i.e. 



Ui = I + X2 - y3, Vi=y2- X3, 

U2 = 1 + Xi - 1/3, V2=yi- X3, 

Us=l + xi-y2, V3=yi-X2, (5.7) 



so that (5.6) simplifies further to 










oo 

>; 

n=— oc 


-p-|- y-'^j)mj+n 


/2/3 
/l 


1111 

i<*<i<3 


{l + x, 

{m- 


•Ji'mi—mj 
3 )mi—mj 


Since (3.28) imph 


.es 












{l + x,- yj)^^_ 

[Vi — ^j)mi-Tn 


-rrij _ 

j 


{l + x 

(%•- 


j yi) ruj —mi 
~ -^i/mj-nii 



{51 



(5.9) 



the symmetry in (5.8) is even larger than is manifested there, namely the full 
permutation symmetry group S3. 



5.2 The two conditions to be satisfied 



The integrable chiral Potts model solution (2.4) of the star-triangle equation 
(2.1) has not six free parameters, but only four independent ones, namely the 
rapidity variables Xp, Xq, x^, and the modulus parameter k. Therefore, two 
conditions must be imposed on the new variables Xj and yj with j = 1, 2, 3. 

The first relation is a linear relation, which is a direct consequence of (4.12) 
and (5.3). It reads 

Xi+X2 + X3 + 2 = yi+y2 + y3- (5.10) 



This is in fact the Saalschiitz condition, which plays such an important role 
in the theory of hypergeometric functions [19,22]. 



25 



The second relation can be found from (4.19), which is derived in Appendix 
A and holds for all A^. From (4.19) we can derive 

AA— A— A— AA A A — A A A — A i^^^\ 

^^pr^^qr ^^pq qp ^^qr^^pr? ^^pq^^qr ^^pr; ^^pq^^qp -■-• \"' -'--'-/ 



From the first equality in (5.11) and (5.3), we then find 

sinTTXi sin7ra;2 sin 71x3 = sin7ri/i sin7rt/2 sin 771/3, (5-12) 



and we may call this the "periodicity condition" due to its relation with the 
periodicity mod A^ property. 

Since (5.12) is a nonlinear relation, we may ask ourselves what the ambiguity 
is in solving X3 and 1/3 from them. Let us use the abbreviations 

sinvrxi sinvrxs , . 

S = - : , T = xi+X2-yi-y2, (5.13) 

smTTi/i sm7rt/2 



which are single- valued functions of xi, X2, yi, 1/2- We then must solve 

ys = X3 + T + 2, sin 7r?/3 = S sin ttx^. (5-14) 

This has the solution 

^3 = ^log g_ei^T ' y3 = x^ + T + 2, (5.15) 



so that the only ambiguity is a translation X3 1-^ 3:3 + M, 1/3 ^-^ ys + M, shifting 
x^ and ys by a common integer M. We will use this freedom below and we 
conclude that we have indeed found the required two conditions for (5.8) to 
hold. 



5.3 A double-sided hypergeometric identity 



Equation (5.8) can be further simplified after we work out j^f'ilfx- We can do 
this using (5.5) and Vj — -Uj = 1 + Xj — yj, for j = 1, 2, 3, which follows from 
(5.7). The result is 

hh ^ TT^ sin TTUi J^ T{yj)T{uj) 

/i sin TTX2 sin ttx^ sin vr (1/1 — Xi) M T{xj)T{vj)T{yj — Xj) 
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Here the functional equation of the Gamma function r{z)r{l — z) = n/ sin nz 
has been used. Substituting (5.16) and (5.7) in (5.8), we obtain 

^ir'oo M r(%- + m, + n) nil nj=i Tiy^ - xj + m^ - m.Y 



where the function G{xi, X2, x^lyi, y2, ys) is given by 

G{xi,X2,X3\yi,y2,y3) = - : z=^ : t r, (5.1^ 



and can also be expressed as a product of ten Gamma functions. 

We can simplify (5.17) further by absorbing the integers rrij in the Xj and the 
yj, i.e. Xj + rrij h-» Xj, yj + ttij i-^ yj. The function G{xi,X2,X3\yi,y2,y3) is 
invariant under this translation. We also note that in this way we utilize the 
freedom in solving the two conditions (5.10) and (5.12) in subsection 5.2. 

As the main conclusion of this section, we find from the large- iV limit of the 
chiral Potts model that the following double-sided hypergeometric identity 
holds: 

^ ^ T{xi + n) _ G{xi,X2,X3\yi,y2,y3) ,^^f.. 



where 



G{xi,X2,X3\yi,y2,y3) = 

3 3 

n r(x,)r(i - X,) n ^{y^ - xi)t{i -y, + x,), (5.20) 

j=2 1=1 

or equivalently (5.18). The identity (5.19) holds provided both the Saalschiitz 
condition and the periodicity condition of subsection 5.2 hold, i.e. 



Xi+X2 + X3 + 2 = yi+y2 + ys, 

sinvra;! sin7rx2 siuTrxs = sinvryi sm7iy2 sinTrya. (5.21) 

Equation (5.19) is clearly also an identity for 



3H3 



a;i,X2,a;3 

1/1,1/2,1/3 



nfTr E n^P-^' (5.22) 



1=1 
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which involves the double-sided hypergeometric function 3H3 as defined e.g. 
by Slater [22]. 

We should emphasize that (5.19) holds whenever (5.21) is satisfied and none 
of the Xj is an integer. The summand and the right-hand side of (5.19) are 
meromorphic functions of their six variables and the sum converges absolutely 
like |n|~^ due to (5.21), see also subsection 4.3. From (5.15) we see that all 
solutions of (5.21) are connected by analytic continuation, and the ambiguities 
in solving (5.21) relate to the way we go around the logarithmic branchpoints 
in (5.15). 

More generally, the N ^ 00 limit treated in this paper corresponds to a 
q = uj ^ 1 limit of the cyclic (basic) hypergeometric functions of [7,24]. 
Looking back, the theory of these cyclic hypergeometric functions seems to be 
nearly synonymous with the theory of the integrable chiral Potts model. The 
fact that we have found yet another new hypergeometric identity confirms this 
point of view. 

No direct proof of (5.19) has been given here. Such a proof should exist and 
the various symmetries of (5.19), 

• Permutations of xi,X2,X3, 

• Permutations of yi,y2,y3, 

• Refiections Xj \-^ 1 — yj, yj \-^ 1 — Xj, for j = 1, 2, 3 simultaneously, 

• Translations Xj h^ Xj + M, yj 1-^ yj + M, shifting Xj and yj by a common 
integer M, for j = 1, 2 or 3, 

• Shifts yi^^ yi + M, yj 1— > yj — M, i j^ j, M an integer, 

should be helpful in such a proof. The left-hand side and the double product in 
the denominator of the right-hand side of (5.19) are both separately invariant 
under the first three of these symmetries, where for the refiection symmetry it 
has been assumed that the conditions (5.21) hold and the functional equation 
of the Gamma function has to be used. The fourth symmetry extends (5.19) 
to a star-triangle equation, as seen above in (5.17), transforming one solution 
of (5.21) into another one. Similarly, the fifth symmetry also maintains the 
validity of (5.19). 

Therefore, the function G{xi,X2,X3\yi,y2,y3) given by (5.18) or (5.20) should 
also exhibit the same symmetries, i.e. 



G{xi,X2,X3\yi,y2,y3) = ^(l-yi, 1-^/2, l-^/sll-a^i, l-a;2, l-ajg) 
= G{x2,xi,xs\yi,y2,y3) = G{x3,X2,xi\yi,y2,y3) 
= G{xi,X2,X3\y2,yi,y3) = G{xi,X2,X3\y3,y2,yi). (5.23) 

This is verified in Appendix D, to which we refer for some details. 
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5-4 Relation with the Dougall-Ramanujan formula 



It is probably now a good idea to ask ourselves how equation (5.19) fits in the 
existing mathematical literature. There are not many identities available for 
double-sided hypergeometric functionsf™] and one obvious source for them is 
the textbook of Slater [22]. It is easily seen that the double-sided hypergeo- 
metric sum pHp({a}; {6}; z) converges only on parts of the unit circle in the 
z-plane, and there are a few identities for z = 1 available. 

Slater gives a general aHa-identity following from a more general sHs-identity 
of Dougall. The identity reads 



3x13 

Ll + a — Xi,l + a — X2,l + a — X3 

r(l — |a)r(l + \a)T{l + |a — xi — X2 — X3) 
~ r(l-a)r(l + a) 

X rr T{l-Xj)T{l + a-Xj) ^ ^^24) 

j=i r(l + a- xi- X2- X2. + Xj)T{l + \a- Xj) 

see (6.1.2.6) of [22], with b = Xi, c = X2, d = X3. Note that this has four free 
parameters, just like in (5.19) where there are two relations (5.21) among the 
six parameters. It is easily checked that (5.19) and (5.24) are different, and 
one may say that (5.24) also has six parameters Xj and i/j = 1 + a — Xj, but 
with two linear relations, whereas one of the relations in (5.21) is nonlinear. 

When a = 0, (5.24) reduces to the Dougall-Ramanujan formula, 

^ ' r{x, + n) _ r{i-x,-x2-x,)ulir{x,) 
ni^jAr{i-x. + n) nn r(i-x.-x,) ' ^- > 

i<«<i<3 



with xi = —X, X2 = —y, X3 = —z in [23] and vahd for ^{xi + X2 + X3) < 1. 
Note, while comparing with (5.19), that we must set yi = 1 — Xi and the 
nonlinear periodicity condition (5.12) is automatically satisfied. Imposing the 
Saalschiitz condition (5.10) implies 

3 3 

Z!^i=2' Y.yi = h yi = l-Xi, (5.26) 
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We are grateful to Professor G.E. Andrews for his comments on this section. 
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so that (5.25) reduces to 

1 r(i-x, + n) ~ ^ '2^ ^y r(xi + i)' 



_£np^?^^-r(i)n^|^. (5,27) 



It is a simple exercise to show that the special case xi + X2 + a;3 = |, ?/« = 1 — Xj 
of (5.19) coincides with (5.27).[3 

Some other interesting hypergeometric sums have been evaluated in [25] , now 
under certain nonlinear conditions, but (5.19) appears still to be new. 



6 Discussion 



There are several reasons why the above results may be of interest, although 
much of this will have to be deferred to future publications. 

First, the integrable chiral Potts model for finite N is intimately related with 
integrable deformations of a series of parafermionic conformal field theories 
[26]. It was originally proposed in [26] that the Fateev-Zamolodchikov model 
[9] constitutes a series of critical lattice models in the same universality classes 
for each A^. This has been numerically checked by Alcaraz [27] for A^ < 8. The 
more general chiral Potts model provides therefore lattice deformations of the 
conformal theory in the chiral-field direction, see e.g. [28,29]. 

This picture should extend also to the large- A^ limit. In section 3 we have 
constructed three different N —* oo limits of the integrable chiral Potts model. 
These lead to three solutions of the star-triangle equations, with an infinite 
sum or finite or infinite integral, as is shown in section 4. Such exact solutions 
with an infinite state space per spin ought to be of interest as they should 
relate to new chiral integrable deformations of parafermionic conformal field 
theories. 

There are not many nontrivial exact results for nearest-neighbor systems with 
infinite spin dimensionality. However, for the chiral Potts model with finite 
A^, exact results exist for the free energy [17,30-32], order parameter [33], 
groundstate energy and excitation spectra [34-36] of the associated quantum 
chain, and surface tensions [37-40]. Many of these results have been obtained 
using systems of functional equations for transfer matrices [41,42]. 

Large- A^ limits of these quantities can be constructed. For example, we can use 



^^ Only the functional relation r(z)r(l — z) = vr/sinvrz and the duplication formula 
T{2z) = 22^-ir(2;)r(z + ^)/r(i) of the Gamma function [19,20] are required. 
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a conjecture for the order parameters [33], i.e. for each n with 1 < n < A^ — 1 
an order parameter is given by the expectation value 

K-) = (1 - A:'^)"^'^""^''''\ (6.1) 



where oq = 0, ■ ■ ■ , A^ — 1 is the random value of a given bulk spin, say at the 
origin. In the limiting regime I, where this spin oq can run through all positive 
and negative integers, (6.1) tends to 



where x is a real number < a; < 27r resulting from the limit 2Tm/N -^ x. 
Similar limits can be constructed from existing exact results for some other 
thermodynamic quantities. 

The precise status of (6.2) must still be determined as the above construction 
presumes the interchange of the N ^ oo limit with the thermodynamic limit 
followed by the field limit defining the order parameter. Furthermore, it would 
be interesting to study all these thermodynamic quantities further within a 
larger oo-state model containing the integrable manifold and we hope to return 
to this in the future. 

Finally, from a mathematical point of view, the N ^ oo limit corresponds 
to a q = uj ^ 1 limit of the cyclic (basic) hypergeometric functions of [7,24], 
which are intimately related with the integrable chiral Potts model. However, 
the connection of the double-sided series (3.26) and (5.19) with more general 
cyclic hypergeometric series will be treated elsewhere. 
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A The Constant A 



In this appendix we verify that the formula for constant A given in (3.2) works 
out in all four cases. 

Substituting (3.3) and (3.4) we get, using (2.17) and (2.22), 
2 cos^|(^p — 0q) 1 + cos ^p cos 0q + sin 6'p sin 0q 

^"^ COS^ I (^q ~ 0p) 1 + COS ^q COS 0p + siu 9^ siu 0p 

^ l + A;^ + 2A;cosgp ^ sin Op sin 0q ^ g27p-27q /^ ^^ 

1 + k'^ + 2k cos 9q sin0psin^q 



and 



sm 



7;i9a — Oo) 1 — COS 9„ COS 6„ — siu ^n sin 6a 



/T \2 2V 1 P' L,ua i/p i^vjo i/q oiii i/p oiii i/q 

^*^ Sin^ I (0q — 0p) 1 — COS 0q COS 0p — siu 0q siu 0p 



(1 + A;2 + 2A; cos ^p) (1 + A;^ + 2k cos ^q 
l-F 

SingpSingq ^ ^27p+27,^ 

sin 0p sin 0q 



(A.2) 



so that we receive agreement with (2.25) and (2.26). The signs of the square 
roots are easily verified for the limit ^q -^ 9p and </)q -^ 0p in case of (A.l), or 
0p ~^ ^p and 0q -^ 9q for case (A.2). 

Similarly, substituting (3.7) and (3.8) into the constants A of (3.2) and using 
(2.21) we have 



^(f) ^ sini(gp-gj ^ ^,,(,^^,^) e^'.-^'^^. - e^p 



P'' sin I (0p - 0q) eT<i+i-^<i - eTp+i-^p 

i(0p+0q-ep-eq)/2 



gi7r(Aq-Ap) 



(A.3) 



A (I) 


COsi(0p-^q) 


hq- 


COsK^ip - 0q) 

-pi(0q-eq-0p + ep)/2 



g7p— i-TrAp g7q+i7rAq i]^ 



(A.4) 



in agreement with (2.27) and (2.28). 
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B Mathematical Details of Large- iV Limit 

In this appendix we prove that 



c.N_>^, sm[7r{j + a-l)/N] 



is asymptotically for large A^ given by 

Sn{a) ^ t^(^y~' [<P^'\^n/N) - 0(^-)(O)] , (B.2) 

where 

0W=/d/log^ (B.3) 



and (p^^^ (z) is its jth derivative. 

The expression (B.l), valid for < n < A^, must be replaced by 

Qf \ V^i sin[7r(j-a)/A^] 

Sn{a) = -2^\og —-— = -SN-n{l-a), (B.4) 

j=i 7r(j-a)/A' 

for —N < n < 0, which has the identical asymptotic expansion (B.2). This 
can be easily verified using 

B,(l-a) = (-l)^B,(a) (B.5) 

and 

0a)(_^) _ 0a)(o) = (-i)^'+i[0(j)(z) - 0(-'')(o)]. (B.6) 

From elementary calculus we note that the derivatives 

d(f){z) sinz d^(f){z) 1 

-- — = log , =cot2:--, B.7 

az z az'^ z 

have Taylor expansions involving the Bernoulli numbers [19,20,43]. We have 
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l^) - Z^ ;/; , i\i - Z^ 



,tt A;(2A;+l)7r2fc' ^""-^^ 

convergent for |2;| < n. Here C(s) is the Riemann zeta function. 
Therefore, we can write 

Sn{a) = j:j:[^) ^U+c^-^Y- (B.9) 

For 

-N <a,n + a-l < N, (B.IO) 



this is absolutely convergent, allowing us to perform the sum over / in terms 
of Bernoulli polynomials [19,20], 

B„(x) = "(x + B)"" ^ J2 (fj BkX--', (B.ll) 



with the properties 

B„(x + 1) =B„(a;)+nx"-\ (B.12) 

Bi+i(«+n) = "(a+n+By+^" = f] ('^+^'| B,(a) x'+l-^ (B.14) 



We find therefore 



^ ' /2vriy B, B,(«) / d^ x'+i 



Z=2 fc=0 



(B.15) 



Here the sum over / converges absolutely within range (B.IO). 

However, the double sum converges only relatively, as one can also verify 
numerically. In order to obtain the result (B.2) we need to interchange the 
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two sums, leading to an asymptotic expansion as a consequence. We find 



^ , , ^ Bkia) N ( d 



A;=0 



k\ n 



dx) frl /•(/ + !)! \n) 

l>2 



[B.16) 



from wliicli (B.2) immediately follows, or equivalently 



N 



Snia) ^ — (j){TTn/N) + {a- i)0^'^(7rn/iV) 



TT 



°° B,+ifa)/7r^' 






0('+i)(7rn/iV)-0('+i)(O)' 



after using Bo(a) = 1, Bi(a) = a - |, 0(0) = ^(^^(O) = 0. 



[B.17) 



The first line of (B.17) is sufficient for our purposes, as we need to keep only 
terms of order 0(1) as A^ ^ oo. Therefore, we conclude this appendix with 
estimating the second line. For n>2,0<a;<l, the Bernoulli polynomial 
can be expressed as [19,20] 



2-n! ^ cos{2rrkx — ^TTn) 
Bn(a;) = — - — — > — 



k' 



{B.li 



so that 



2-n' 
Bn{x)\<——an), forn>2, < a; < 1, 



(B.19) 



with 1 < ({n) < C(2) = 7r^/6. Also, for m = 2, 3, ■ ■ ■ , and \z\ < n, we have 
from (B.8) 






2fc+l-m 



[B.20) 



where [x] is the smallest integer > x. Since all the coefficients in (B.20) are 
negative, we have 

0(™)(z) - 0(™)(O) = (-l)'"+^[0(™)(-z) - 0("HO)] < (B.21) 



and monotonically decreasing for < z < n. Therefore, we can estimate 
(B.20) by replacing all Ci'^k) by its minimum 1 or its maximum (^(2[|m]), 
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which makes the sum a binomial-type expansion. We find 

^M(^) < _(sign2)™+^ U('")(2) -0(™)(O)j < C{2\lm])x^"'\z), (B.22) 



where 



^(2n+i)(^)^(2n-l)! 



1 



1 



(7r-|2|)2«-i (7r+|2|)2"-i 
1 1 



{n-\z\y^ 
which are both positive for n = 1,2,- ■ ■ 



U+lzlV'^ 7r2" 



(B.23) 



Bounds (B.19) and (B.22) are quite sharp, as can also be numerically verified. 
The bounds show that (B.2) is an asymptotic expansion, with absolute value of 
terms roughly bounded by 1\/{'kNY and error less than the bound on the first 
ignored term. This is true for N >2 and \n\ < N. But the expansion becomes 
particularly useful for large A^ and the full range of n taken as \n\ < ^N. 



C Derivation of (4.8) 



The factor F = Fpq is obtained substituting (4.4)-(4.8) into (4.3). We can 
simplify it using (4.1) and the functional equation of the Gamma function 
r(a;)r(l — x) = 7t/ sinirx. This gives 



F-2 f—Y'^ A^I'^A^^^''' r(a)r(/j)r(l-/j + a)sin7raffl 

V27rj AV2 r(«)r(/?) 



(C.l) 



Part of this expression can be simplified further using A = sinTr/^/sinTTO;, 
(4.8) and (4.1). More precisely, 



4A A^^'> sin^ Tra'^'^^ = 4A sin Tra^''^ sin tt/S^'^^ 



(f)^ 



2 A cos 7r(a(^) - ^^^^) - 2 A co^ixi^^^ + /3 
-2A cos 7r(a - ^) + 2A cosh log A 

sin^7r(/9 — a) 



A^ -2Acos7r(a -/3) + 1 



sin^ Tia 



(C.2) 
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Taking the square root of this in domain (4.2) and applying it to (CI) we find 

^ / AT N /^- r(a)r(;^)r(l-:^ + a)sin7r(^-a) 

\2nJ AV2r(a)r(/5)sin7ra ' ^ ' 



This can be further reduced to (4.8) using the functional equation of the 
Gamma function again. 



D Verification of (5.23) 

In this appendix we show that under the conditions (5.21), or 

3 3 3 3 

'^ + J2^j = J2 Vv n ^^"^ ^^j = n ^^"^ ^yj = ^/(2i)^- (D.i) 

j=l j=l j=l j=l 

the symmetry relations (5.23) hold or, equivalently, 

sin7r(yi -Xi)sin7r(i/2 - a;i) sin7r(y3 - Xi] 
sin nxi 



a/{2\f (D.2) 



is fully symmetric both in {xi, X2, x^} and in {t/i, 1/2, 1/3} and it is also invariant 
under Xj 1-^ 1 — y^, yj \-^ 1 — Xj, for j = 1,2,3 simultaneously. Indeed this 
gives (5.23), as G{xi,X2,X3\yi, 1/2,113) = (27ri)V/(T. 

With the definitions 

^i = e'"^% 7]i = e'"^% for i = 1, 2, 3, (D.3) 

the conditions (D.l) become 

3 3 

^1^2^3 = VIV2V3 ^ P, Ili^^-^7') = IliV^-V^')=r. (D.4) 

i=l i=l 

We can now expand pr and rearrange terms. We obtain 

P^E(r-^r^) = E(c^-^n. (D.5) 

j=l 4 = 1 

We use this to expand a. Successively, we find 
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P(e? - 1) 

= E(^--C')/P- (D.6) 

i=l 

This shows that a has the required permutation symmetries. The invariance 
of a under the reflection symmetry 

6 ^ -^»"\ Vi ^ -C^ (^ = 1. 2, 3), and p h-> -p-\ (D.7) 

follows from (D.5) and (D.6). As r obeys these symmetries trivially, we can 
now complete the proof of (5.23). 
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